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We discuss a supergravity description of the metastable state that is created by 
a stack of D3-branes placed at the tip of the KS background. When the number p 
of the D3-branes is large g s p ^> 1 the characteristic curvature of the corresponding 
gravity dual is large in stringy units and one may expect the background to be regular 
everywhere. Starting from the distances of order R ~ (gspY^a' 1 ^ 2 away from the tip 
the new background can be well approximated by a linear perturbation around KS. 
By applying the appropriate boundary conditions in both IR and UV we found the 
lowest KK mode of the corresponding linear perturbation. The solution we found 
contains VEVs of the SU(2) x SU(2) invariant operators at the linear order in p. 
As a non-trivial check we calculate the ADM mass which exactly matches the probe 
approximation. As a byproduct we also found a gravity background dual to the KS 
theory deformed by the operators W 2 and W 2 W 2 with small coefficients. 
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1 Introduction and Summary 

The SUSY-breaking metastable states [1] is an interesting topic which plays im- 
portant role in different phenomenological applications. One famous example is the 
metastable state in SU(M(k + 1) — p) x SU(Mk — p) theory found by Kachru, Pear- 
son, and Verlinde 0. They start with the SU(M(k + 1)) x SU(Mk) theory dual to 
the deformed conifold solution of Klebanov and Strassler |3j and place p D3-branes 
at the tip of the conifold. As a result of non-abelian dynamics the stack of D3-branes 
blow into a NS5-brane. The NS5-brane carries p units of D3-brane charge. The 
corresponding SUSY-breaking state turns out to be metastable when p/M is smaller 
than a certain critical value of order 0.08. This state can tunnel into the true vac- 
uum with NS5-brane disappearing and M — p D3-branes emerging. Although the 
field theory description of the metastable vacuum is out of control due to strong cou- 
pling, the schematic picture is the following. The original SU(M(k + 1)) x SU(Mk) 
theory has a rich space of vacua which includes the baryonic branch, described by 
the KS solution, and the mesonic branches labeled by I > 1 described by the KS 
solution with IM mobile D3-branes placed on the conifold. After placing additional 
p D3-branes the mesonic branches survive while the baryonic branch gets uplifted. 
Thus the metastable state of the SU(M(k + 1) — p) x SU(Mk — p) theory is the 
former baryonic branch of the SU(M(k + 1)) x SU(Mk) uplifted in the presence of 
the D3-branes. 

Since the field theory does not provide a reliable description of the metastable 
state it is desirable to describe the corresponding physics with help of the gravity 
dual background. A priori we do not know if supergravity approximation will be 
reliable near the location of D3-branes, but in any case sufficiently far from the tip 
the supergravity description should be reliable. Therefore, no matter how complicated 
the dynamics is in the direct vicinity of the D3-branes, in the end it should boil down 
to some boundary conditions for the supergravity fields at the location of D3-branes. 
These boundary conditions come from the probe action for the D3-branes for p ~ 1 
or NS5-brane for p ^> 1. To illustrate why the probe action leads to the correct 
boundary conditions despite possible complicated dynamics we turn to the classical 
electrodynamics. The point charge q creates a potential Aq = q/r that is singular 
at the origin. The classical electrodynamics breaks down at the distances shorter 
than the inverse electron mass m" 1 . One has to rely on QED beyond that scale. Yet 
classical electrodynamics is perfectly valid for the distances much larger than m" 1 . 
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Moreover to find the correct boundary conditions for A at the origin one does not 
have to know QED. Rather the boundary conditions come from the probe action 
J d 4 xA fl j^ l that can be established by probing physics at small energies i.e. distances 
larger than m" 1 away from the charge. The story with the D3-branes is no different. 

Although the task of finding the supergravity solution dual to the metastable 
state is well-defined, technically it is very challenging. First few steps were taken 
in the original paper [2] and later in [3J. The qualitative picture is the following. 
Very close to the D3-branes, at the scales much smaller than the curvature scale of 
the KS solution R 2 KS ~ (g s M) l l 2 a' the geometry looks like an AdS$ throat of radius 
~ (g s p) l l 2 a' created by the D3-branes. If p ~ 1 the corresponding geometry 
is highly curved and no supergravity description valid everywhere is available. But 
when g s M ^> g s p ^> 1 the curvature radius is large in the stringy units and the 
solution is expected to be regular everywhere. The AdS$ throat created by the D3- 
branes is corrected by small perturbations coming from the KS backround. Some of 
these perturbations break (anti)SUSY of the D3-branes. There are also modes, such 
as the ISD flux, which are relevant in the dual field theory sense. The relevant modes 
increase as we come closer to NS5-brane and may, though not necessarily, destroy the 
geometry. If we focus on the effect of the ISD flux and neglect the effect of other modes 
the deep IR region near the NS5-brane will look similar to the Polchinski-Strassler 
solution [5]. In fact the latter is exactly the IR region of the AdS§ created by the NS5- 
brane carrying the D3-charge, perturbed by the relevant AISD mode. Hence we expect 
the full supergravity background dual to the metastable state to schematically look 
as follows. In the very IR near the NS5-brane the solution will look somewhat similar 
to the (anti)Polchinski-Strassler solution. As we move away from the NS5-brane the 
Polchinski-Strassler solution turns into the AdS$ throat perturbed by some relevant 
modes with small coefficients. The AdS$ is glued to the tip of the Klebanov-Strassler 
solution forming a smooth weakly curved background. Clearly to find such a solution 
beyond our reach at this point. Nevertheless we may not need the full solution to 
describe interesting physics. If we are interested in the vevs of various operators in the 
metastable state or in numerous phenomenological applications of SUSY-breaking, it 
is enough to know the behavior of the solution at large radii r 3> -Rd 3 . I n this regime 
the solution in question is only slightly different from the original KS background and 
can be found in linear expansion around the latter. In the linear approximation we 
will loose quadratic and higher in p effects e.g. the difference between the stack of D3- 
branes and the NS5-brane. Yet this approximation contains all leading in p effects e.g. 
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SUSY breaking. Moreover working at the linear level admits further simplifications: 
it is self-consistent to focus only on the SU(2) x SU (2) invariant modes. Let us point 
out that such a solution is a part of a more complicated linearized solution produced 
by the localized D3-branes at the tip, not a solution produced by the D3-branes 
smeared around the tip. The smeared configuration is not stable as the D3-branes 
attract each other. 

Up to date there were several works attempting to construct the SU (2) x SU(2) 
invariant part of the linearized solution. The pioneering work [6] used the expansion 
around the KT solution [7j to identify the modes produced by the D3-branes. Since 
the KT solution is singular at the origin, the proper boundary conditions in the IR 
could not be imposed and the resulting solution was valid only up to the leading 
order in 1/r 4 . Although this accuracy is enough for many interesting applications 
it is not enough, for example, to determine the force induced on a probe D3-brane 
placed in the background. Another work [8] attempted to find a solution expanding 
in Taylor series near the tip. Such a solution can not be continued to the UV region. 
As a result half of boundary conditions can not be imposed and the resulting solution 
can not be identified. The omission of the approaches based on the perturbative 
expansion in small or large radius was intended to be cured by the approach taken 
by Bena, Grana, and Halmagyi in [H] (see also [TD] and [TT] for further development). 
They used the technique developed by Borokhov and Gubser in [T2] and latter refined 
in [13] . The idea is to make a change of variables representing seven second order 
equations governing linear perturbations around the KS background through fourteen 
first order equations. This allows to find the solutions valid at any radius, albeit in 
an indirect form. The authors of [H] found that the boundary conditions for the stack 
of D3-branes will necessity require the solution to be singular at the tip. While the 
singularities of the warp factor and the RR 4-form directly follow from the boundary 
conditions provided by the D3-branes and were expected, the singularity of the 3-form 
flux was confusing. Indeed the seven modes system describing linearized deformations 
of the KS background has fourteen solutions: seven should be regular in the IR while 
singular in the UV and seven should be regular in the UV while regular in the IR. 
This logic suggests that one should be able to pick any boundary conditions in the 
IR at will, at least if the UV behavior is not restricted. Yet it was emphasized in [9] 
that once the solution contains the singular mode directly related to the presence of 
the D3-branes the 3-form flux becomes singular as well. The same observation was 
also made in [8] . Since the D3-branes directly source only one IR singular mode (the 
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combination of the warp factor and the RR 4-form $_ in the notations of [13]) that 
does not include any modes related to the 3-form flux, why then one can not choose 
the 3-form flux to be regular? This is due to the coupling of $_ to the flux which is 
evident from the following equation of motion [T3] 

dr 

(d + i- — ARe)($_G+ + $+G'_) = . (1.1) 
Imr 

Here G± are the ISD and IASD parts of the 3-flux while $+ is yet another combina- 
tion of the warp-factor and the RR 4-form. $ + couples to the D3-branes and stays 
regular in their absence. In the KS case $_ = G_ = while $ + , G + ^ 0. After intro- 
ducing the D3-branes, $_ calculated at linear order diverges at the origin. Hence to 
compensate this divergence in (11.11) there should be a singular flux G-. Let us note 
here that at the non-linear level very close to the D3-branes $_ will be proportional 
to the inverse warp-factor $_ ~ h~ l and hence regular. That is why the reason for 
the IASD flux to be singular will disappear in a full non-linear solution. 

The aforementioned singular mode of the 3-form flux is in fact an IR regular mode 
of the linearized system |15| . In particular it is normalizable near the tip. There is 
another mode of the 3-form flux which is more singular at the tip but is UV regular. 
This mode may not and should not be turned on in the solution describing the D3- 
branes placed at the tip. Here we outlined the behavior of the solution near the tip 
only schematically. Later in the text we will return to this point. 

As a side note we would like to emphasize that the singularity of the 3-form flux 
discussed above is not an effect of smearing. The SU(2) x SU (2) invariant linearized 
solution is the lowest KK mode of the linearized solution sourced by the localized 
D3-branes. Therefore if the former is singular so is the latter. This is also evident 
from (II. ip as the linearized is singular in the localized well. Rather the 

singularity of the 3-form flux is an artifact of linearization. In any way we can only 
trust linear solution to describe physics for the distances r 3> (^p) 1 / 4 ^' 1 ^ 2 away from 
the tip, where the solution is regular. 

In this note we construct the SU(2) x 377(2) invariant mode of the linearized 
solution describing the metastable state of SU(M(k + 1) — p) x SU(Mk — p) theory 
in the leading in p order. We further develop the approach of [5] and show that there 
is a unique solution that satisfies the proper boundary conditions in the IR and the 
UV. We found this solution expressed indirectly in terms of some integrals. We use 
this solution to calculate the ADM mass of the metastable state. It turns out to 
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be equal to the probe value 2pT 3 h~^ s (0). As the ADM mass is sensitive to the IR 
boundary conditions this provides a non-trivial check of our result. As a by-product 
we also obtain a gravity dual of the KS theory perturbed by the operators W a W a and 
W a W a W l3 Wi3 with small coefficients. The resulting solutions are regular everywhere 
and can be used to study gauge theories with softly broken SUSY. 

This paper is organized as follows. In Section 2 we review the system of linear 
equations governing the perturbations around the KS background and find a general 
solution parametrized by 14 integration constants. In Section 3 we formulate the 
boundary conditions for D3 and D3-branes and construct the SU{2) x SU (2) invariant 
mode of the linearized solution describing the metastable state. In Section 4 we 
calculate the ADM mass of the metastable state. In Section 5 we briefly describe the 
gravity background dual to the KS theory perturbed at linear level by the SU(2) x 
SU(2) invariant, parity and X-symmetry even operators. 



2 System of Linearized Equations 



We start with a system of coupled ODEs introduced in [T2] and follow the nota- 
tions of [9]. The system is specified by the kinetic term 

G ab d<j) a d(f) b = 

e 4p+4A (dx 2 + hy 2 + 6dp 2 - 6dA 2 + ± e -*- 2x (e~ 2y df + e^dk 2 + 2e 2q, dF 2 ) + \d§ 
and the superpotnetial 

W((f>) = e AA - 2p ~ 2x + e AA+4p cosh y + I e 4A+4p ~ 2x (/ (2P - F) + kF) . (2.2) 

The solution can be expanded = (p + 5<p + .. where O solves the unperturbed 
system 



b a Q G ab dW 



dt 2 d(j) h 

The equations for 5<j) could be linearized 

d£a 



(2.3) 

4>=4>o 



(Jt - &M b a , (2.4) 



- G a % , (2.5) 



dt 
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where 

1 9 / dW 



(2.6) 



The unperturbed solution describing the KS background (we put g s and the de- 
formation parameter e to be equal to one) is 



e 



x ° = ^ /i 1/2 (cosh t sinh t-t) 1/3 , 



e 



e yo = tanh(t/2) , 

6po 24 (cosh t sinh t-t) 1 / 3 
h sinh 2 1 



e 6A ° = — ^— /i(cosht sinh t - t) 2/3 sinh 2 1 , (2.7) 
3 • 2 9 

(tcotht- l)(cosht- 1) 



fo = ~P 



sinht 

(fcothf- l)(cosht + 1) 



Fo = P 



sinht 
(sinh t — t) 



sinht 
$o = 0. 

Here P = —\Mo l and the warp-factor is defined as 

^ _ e -4A -4p +2a;o 

= h - 32P 2 [ T tCOtht ~ 1 (cosh t sinh t-tf'^dt. 
Jo sinh t 

The constant h = h(0) is chosen such that h approaches zero at infinity. 

To write down the linearized equations for <50, £ we change the variables 0, £ — >■ 0, £ 

= (x - 2p - 5A, y, a; + 3p, x - 2p - 2A, f, k, F, $) , (2.8) 

1 = (36 - 6 + 6, 6, -36 + 26 - 6, -36 + 6 - 26, 6 + 6, 6 - 6, 6, 6) , 

while the inverse transformation is 

= (x,y,p,A,f,k,F,*) (2.9) 

= (^(-201 + 203 + 50 4 ), 02, ^(201 + 303 - 50 4 ), ^(-01 + 04), 05, 06, 07, 0s) , 

e = (|(36 + 6 + 6), 6, 6 +6, -6 -6, |(6 + 6), |(6 -6), 6, Is) . 



The linearized equations of motion are (from now on we denote the perturbation 



part 5(f) simply by 0) 

g = e- 2xo (2Pf -F (f -k ))i 1 (2.10) 

g = -e^U2P/o - Wo - MKi (2.11) 

g = (2.12) 

g = -It- V 2 *»(P-P )g (2.13) 

g = -sinh(2yo)|5-cosh(2y )|6 + ^e- 2a:o (/o-A;o)| 1 (2.14) 

g = (Pe 2 ^-sinh(2y )Po)g + (Pe 2 ^- cosh(2y )P )g + |(/o - fc )# 7 (2.15) 

g = 3e- 2 ^- 6 ^3+(5e- 2 ^- 6p «-e- 2 ^(2P/ -F (/ -A; ))6 (2.16) 
g = g cosh yo + ^ sinh y (2g + g + g) 

+2 ((Pe 2 ^ - cosh(2y )Po)g + (Pe 2 ^ - sinh(2t/ )P )g) (2.17) 

and 

0' 8 = -4e- 4(Ao+po) g (2.18) 

0' 2 = -coshy o 2 -2e- 4(Al,+Po) g (2.19) 

3 = _ 3e -6po-2.o0 3 _ sinh yo fc _ l e -4(A 0+P o)( 9 | i + 5 | 3 + 2 g) (2.20) 

0; = 2e- 6po - 2 >3 - sinh + I e - 4 ^° + ^(g + 3g) (2.21) 

o 

0' 5 = e 2w (P - 2P)(20 2 + 8 ) + e 2w 7 - 2 e - 4 ( Ao+po ) +2 ^ ,+ ^(g + g) (2-22) 

0' 6 = e - 2 ^(Po(20 2 -0 8 )-0 7 )-2e- 4 ^ + «') +2 ^-^)(g-g^ (2.23) 

4>'r = \ (05 - 06 + (fco " /o)0 8 ) - 2e^'+«') +2 ^g (2.24) 

01 = V 2a;o (/o(2P-Po) + A;oPo)(20 1 - 203 - 504) + ie- 2;Co (2P-Po)05 
z 

+ ie- 2 ^P o 6 + ^e- 2 *°(fc - / o )0 7 - i e - 4 ^+^)g (2.25) 
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2.1 Solving £ equations 

To solve the equations we follow the basic procedure outlined in [9]. First we solve 
the equations for ^ and then deal with 0j. We notice that the equations (12.101) can 
be integrated 

6 = X x h , (2.26) 

where X\ is some constant. The reason why the expression for £i is so simple will be 
explained below in Section 3. The KS warp-factor h is not an analytical function and 
since £i enters the equations for all other £j we won't be able to solve them explicitly. 
Therefore we will present the solutions for £j assuming Xi — first and restore X x 
later. The solutions parameterized by the constants X 2 , . . . , X 8 are 

& = 0, (2.27) 
£4 = X 4 , 

~ X 6 X 7 (cosh(t)sinh(t) -t) X 5 t 
S6 — ■ , 77; "I rrTT; K 



2sinh(t) sinh(t) sinh(t) ' 

~ X 6 cosh(t) X 7 (sinh(t) cosh(t) 2 + t cosh(t) - 2 sinh(t)) + X 5 (t cosh(t) - sinh(t)) 
sinh(t) sinh(t) 2 sinh(t) 2 

~ _ YMPY X 6 P(tcosh(t)-sinh(t)) 
£ 8 - X 8 + PX 5 — — 

(X 7 - X 5 )P(tcosh(t) - sinh(t))(cosh(t) sinh(t) - t) 



sinh(t) 3 

| 3 = 4X 3 (cosh(t)sinh(i) - 1) , 

~ . 4X 3 tcosh(t) X 4 cosh(t) X 5 P(sinh(t) 2 - 2 sinh(t) cosh(t)t + t 2 ) 
£2 — 2X2 smh(t) 1 1 



sinh(t) 3 

X 6 P(cosh(t) sinh(t) - t) X 7 P(2 sinh(t) cosh(t) 3 t - sinh(t) 2 - t 2 ) 



2sinh(t) 3 sinh(t) J 
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Now, if we restore X 1 the solutions above will shift as follows & — >■ & + X ± Q 
Ci = h, (2.28) 
C 4 = -h , 



5 



6 



where 



16P r da; 



(cosh(t) sinh(t) - tfl 3 ' 

(cosh(t) Sinh(t) - t) /' gg| - /* (^H-)^Hx)-x)G(x) dx 



2 sinh(t) 

16P t 

3 sinh(t) (cosh(t) sinh(t) - t) 2 / 3 ' 



P(tcosh(t) - sinh(t)) 
smh(i) 2 

16P /"* / P:r(:rcosh(:r) - sinh(x)) 
ax 



+ 



sinh(x) 3 (cosh(x) sinh(:r) — x) 2 / 3 
P(sinh(a;)(cosh(x) 2 + 1) - 2xcosh(a;)) f x dy 



f 



( 3 = -h+ ^(cosh(t)sinh(t) -t) J 
( 2 = sinh(t) J dx — 



sinh(x) 3 J (cosh(y) sinh(y) — y) 2 ! 3 ) 

2smh.(x) 2 h(x)dx 



(cosh(a;) sinh(a;) — x) 2 
' 2 (cosh(x) sinh(:r) — x) f x 2sinh.(y) 2 h(y)dy 



9 sinh(x) 2 J (cosh(y) sinh(y) — y) 2 

32P 2 t(cosh(a;) 2 + 1) - 2 cosh(a;) sinh(x) f x dy 



+ 



3 sinh(:r) 4 J (cosh(y) sinh(y) — y) 2 / 3 

2P(sinh(a;)(cosh(x) 2 + 1) — 2xcosh(a;)) f x 



smh(x) 4 



/X 
Ce(y)dy 



G ^_ 16P f 4tsinh(t) + 2cosh(t) f f dx 



r dx \ 

J (cosh(t) sinh(t) - t) 2 / 3 J 



3 V3(cosh(t)sinh(t) -t)5/3 S inh(t) 2 J (cosh(t) sinh(t) - t) 2 / 3 

The zero energy condition ia-jj 1 — implies 
d(x + p-A) ~ d(x/3 + p) ~ d(x/3 - A) ^~dy ~ d(f + k) ~ d(f - k) ~ dF 

6 It + 6 — It — + u — It — + ^Tt + ib ^2dT~ + 6 ^T~ + 6 dF 

= ~X 3 + 2X 2 - PX 5 + 3PX 7 = . (2.29) 
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2.2 Solving equations 

Now we can proceed with the equations for 0. Although we will find a general 
solution, we are primarily interested in the solution dual to the metastable state in 
the SU(M(k + 1) — p) x SU(Mk —p) theory. The latter must have the same behavior 
at UV as the original KS solution (besides the shifted D3-charge). That will help us 
to fix seven integration constants by requiring that all UV singular modes are turned 
off. Besides, we will impose regularity at the IR of all modes except the one that is 
directly coupled to the D3-branes. That will fix another six integration constants. 
The boundary condition for the remaining singular mode will be fixed in Section 3. 

The first equation for the perturbation of the dilaton 8 is easy to solve 

r°° 64 

^ = Y 8 + , . fA . UA , W3 e 8 • (2.30) 
J t (cosh(t) smh(t) — tp 6 

We choose the integration constant Y s = such that 0§ — > at infinity i.e. there is 
no source for the operator Tr(F x 2 + P 2 2 ). The convergence at IR implies 

X 8 = ^(X 6 -6X 5 ) . (2.31) 
6 

The equation for <p 2 is also straightforward 

32 [* smh(x)£ 2 (x)dx | Y 2 



sinh(t) J (cosh(x) sinh(x) — x) 2 / 3 sinh(t) 
Clearly to avoid 1/t singularity in the IR we put Y 2 = and also require 

/2\ 1/3 2 

-16PM -J X 1 + -h X l + PX 6 -X 4 = (2.33) 

to avoid \og(t)/t behavior at small t. In general (f) 2 defined by f)2.32p goes to zero as 
e~*/ 3 at large t. Given that the fluctuation 2 is dual to the operator Tr(AiAi + A 2 A 2 ) 
of dimension 3 we must require the leading asymptotic e - */ 3 to vanish. As follows 
from f)2.32p this would require the leading term of order e* in £2 to vanish. As there 
are two terms of the same order, te l and e*, this implies two constraints 

- lx 3 + PX 7 + aX 1= 0, X 2 + + 0X X = , (2.34) 
3 6 

where (at+/3)e* is the UV asymptotic of £2- The coefficients a and (3 can be calculated 
numerically. Their values are not important for what follows. 
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At first glance it looks surprising that the condition that there is no mass for 
Tr(AiAi + A2A2) leads to two, rather than to one condition. In fact besides Tr(AiAi + 
A2A2) there is another operator Tr(AiAi — A2A2), also of dimension 3. Although 
these operators are independent, the corresponding modes mix in the KS case. The 
constraint (I2.34p ensures that both perturbations are absent. 

At the next step we calculate </>3 

v. r f ( cos H%) sinh(x) -x) ~ 

Y * + l { s^w 02 ~ (2 ' 35) 



(cosh(t) sinh(t) - 1) 

- jj(cosh(x) sinh(x) - x) 1/3 (7£i + 5£ 3 + 2X 4 )^ dx 

Again, the integration constant Y 3 = should be chosen to avoid the 1/t 3 singularity. 
To cancel 1/t behavior we must require 

2N l/3 



X 6 - 16P X x = 0. (2.36) 

The \og(t)/t term vanishes automatically at this point due to (12.331) . The leading 
UV asymptotic of 3 ~ e 2 */ 3 is due to e 2t divergence of £ 3 . This mode is dual to the 
perturbation by the bottom component of the operator TrW 2 W 2 of dimension 6 and 
hence must be canceled by an appropriate choice of X 3 

A 3 + 7A1 = . (2.37) 

The coefficient 7 can be calculated numerically and, again, the exact value is not 
important for what follows. This constraint assures the coefficient in front of the 
leading term e 2t in £ 3 vanishes. In turn the first constraint (12.341) implies that the 
coefficient in front of the leading term e* in ^6)^7 vanishes as well. 

The integration of gives 

7 P , ( 4sinh(a;) 2 ~ , , fafx) 

01 = Y 1 - / dx[ - 1 Mx) + 4^ 2.38 

J t \ 3(cosh(x) smh(x) — x) smh(x) 

_ '~3(cosh(x) sinh(x) — x) 2 / 3 ^ 4 ^ 

We choose the integration constant Y\ — such that <pi goes to zero at infinity. In 
fact it is a gauge choice to keep <\>\ vanishing at infinity. Finite value of 0i would 
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only lead to a non-physical rescaling of the warp-factor A. As the IR log(t)/t and 1/t 
behavior vanishes due to (I2.33P and (I2.36P and the log(t) behavior vanishes due to 
the zero energy constraint (12. 29ft . we do not get any new constraints from regularity 
of 0i at small t. 

The equations for 5 , 6 , 07 can be combined in such a way that the dependence 
on 05,06 drops and one gets a second-order differential equation for 7 . The corre- 
sponding solution depends on two integration constants 

= (cosh(t) sinh(t) - t) fr + /o ggg^ggg 

^ 2sinh(t) 

Y 7 ^ F 7 (cosh(t) sinh(t) - t) 



sinh(i) sinh(t) 

= 16P(tcosh(t)-sinh(t)) - 2- 
7 ~" sinh(t) 2 (cosh(t)sinh(t) -t) 2 /3^ 8 3 ^ 

2P(smh(t)(cosh(^ + 1) -2tcosh(t))- 2 _ ^ 



sinh(t) 3 

The choice of Y? — is dictated by the 1/t singularity in the IR. In the case of general 
02 the leading asymptotic of Gj is e - */ 3 and therefore to cancel the e* growing term 
in 07 we required to put Y7 = 0. The corresponding perturbation of dimension 7 is 

jd 2 ew 2 w 2 . 

The integration of 05, 06 is straightforward now 







(2.40) 



As we started from three first order equations for 05,6,7 only three combination of 
Y 5 ,Y e ,Y 7 , Y 7 are independent. There is a constraint that expressed Y G through Y 5 

16 „ /2 N 1/3 



n = Y 5 - —P ^-j h X 1 . (2.41) 

The fact that Y 5 = 05 (0) is not equal to Y 6 = 05 (0) leads to the divergent NS-NS 
field H 2 ~ t~ 2 at small t. As is seen from (12.411) this singularity is always present for 
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Xi 7^ i.e. when the D3-branes are present [9]. In fact regular, although non-equal 
at the origin 5 and 06 corresponds to the IR regular mode of the linear system in 
question [15]. In particular this mode is normalizable in the IR. There is also an IR 
singular mode with 05,6,7 diverging at the origin. That mode should not be turned 
on as there is no direct delta-function source for this mode coming from the probe 
action for the D3-branes. 

As was mentioned above, absence of the perturbation by the operator Tr(AiAi — 
A2A2) is already covered by the constraint (j2.34j) . Indeed in (12.37P we required that 
the leading e 2t asymptotic of £3 vanishes. Therefore the only term potentially leading 
to te t behavior of £2 (and correspondingly to te~ t//3 behavior of 02) can come from the 
e* divergence of £g- But the first constraint in ( I2.34p is nothing but the requirement 
that and hence £7, does not have the leading e* term. As a result the fluctuation 
07 given above goes to zero as e - ', not as e~'/ 3 . This, together with the equation 
(I2.24p implies that 05 — 06 approach zero at infinity also as e~*. At the same time 
05 + 06 at infinity is some constant that is related to the net D3-charge. In general the 
D3-charge is given by — (/(2P — F) + kF). Therefore shifting / + k at infinity 
by a constant while keeping f — k,F the same corresponds to changing the total 
D3-charge resulting from adding Nd3 D3 or A% 3 D3-branes 

f(r -> 00) - f KS (r -> 00) = k(r -> 00) - k KS {r ->• 00) ->• -^Vds - i% 3 ) ( 2 - 42 ) 

The difference f(r —> 00) — fxsij" ~ * °°) is n °t quite the same as 5 (t — > 00) as the 
relation between the radial coordinate r and the dimensionless coordinate t involves 
the deformation parameter of conifold e: r 3 ~ ee*. Since the deformation parameter 
could be different for the original KS background (we choose it to be e = 1) and for 
the background perturbed by the D3-branes the condition (I2.42p will determine e as 
will be discussed in the next section. 

The last remaining mode is 

f°° (1 ~ 16A h 2 

Y<-J (2.43) 



+ 8P 



;sinh(x) + x)0 5 + (sinh(a;) - x)0 6 - 2 xcosh s ^ ) h( g nh(a:) ( 
sinh(x)(cosh(x) sinh(x) — x) 2 / 3 



Clearly the integration constant Y± must be put to zero. Otherwise the leading e 44//3 
behavior of 04 would correspond to the perturbation by the operator J d 2 6d 2 6Ti(W 2 W 2 ) 
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of dimension 8. The IR behavior of 04 is singular 

Sh = 2hil = v (1) 2/3 jWBlzjg&Q) + (2 . 44) 

3 V 2 J £ 

because of the delta-function at the origin associated with the presence of iV D3 D3 
or iVp 3 D3-branes [9]. We can relate the coefficient in front of the singular term in 
f)2.44p to the total number of Ad3 D3 or A% 3 D3-branes using the equations of motion 
which can be found in the next section 

^ (X 1 h(0) 2 - 3PY 5 ) = N D3 + N m . (2.45) 

This constrain can be understood as an equation for Y5. 

Thus we have fixed all integration constants except for X\. We will do it in the 
next section. 



3 Solution for D3-branes 

As a warm-up we start with the solution for the D3-branes placed at the tip of 
the conifold i.e. put N^ 3 = in the boundary condition ( I2.42p . As a result we have 
05, 06 shifted by a constant 

2ttP 



, }p Nm. (3.1) 

At the same time the D3-branes are mutually BPS with the background and therefore 
experience no force. This fixes X\ (and hence all other constants) to be zero [9] . The 
only non-trivial quantity is the perturbation of the warp-factor 



00 



5h = 2/i04 = 4W 2 iV D 3 / . u r d l ( r . (3.2) 

(cosh^x) smh(x) — x) 



All other functions 0j for % 7^ 5, 6, 4 are zero while Y 5 = 05 (i) = 4>§{t) = ~ igr^D3 in 



2-kP 

lui t 7= <j,u,4: aits ztsiu wmie 1 5 — ip^yvj — <±>q\^) — ~ 

agreement with (I2.45p . 

The formalism used in this paper assumes the gauge choice F5 = F3 A B 2 . Hence 
the fields /, k and their perturbations 05^ control not only B 2 but also F§. In fact we 
can make a gauge transformation shifting /, k in B 2 by a constant while keeping F5 
intact. It is convenient to choose a gauge such that the integral of B 2 A F 3 over T 1,1 
vanishes at the origin. Then the value of F 5 at the origin defines the Page charge 

Q™ = 7TTlvl F 5 -B 2 AF 3 . (3.3) 
{A-K z a'Y J T i,i 
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For the small perturbation around the KS solution the Page charge is given by Qdz — 
— -^pY 5 . The Page charge is not gauge invariant. When the integral 

B 2 (3.4) 

s 2 



47TV 

becomes larger than 1 one should make a gauge transformation shifting (13. 4p back 
to zero. This would change the Page charge by M. On the field theory side this 
corresponds to the Seiberg duality which changes ranks of the gauge groups. But the 
Page charge is quantized and mod M conserved and therefore it counts (mod M) the 
number of D3 or D3-branes placed at the origin. 

If we start with the KS background and add Nds = IM, I > 1 D3-branes the 
corresponding solution would be a gravity dual of the Z-th mesonic branch of the 
SU(M(k + 1)) x SU(Mk) theory [16J. To have the same behavior in the UV as the 
original KS solution the new solution must have different deformation parameter e. 
To find the relation between e on different branches we can look at the behavior of 
k, f at infinity. For the original KS solution dual to the baryonic branch of moduli 
space we have 

/, k _> - P (t - 1) = -P(logr 3 - log ejj - 1) , (3.5) 

while for the Z-th mesonic branch 

/; k _ i + ^) = _p(i ogr 3 _ log e 2 _ ! + ^) . (3 . 6) 

The condition that (13. 5p and ( 13. 6 p coincide gives 

e\ = e 2 e^ (3.7) 

in agreement with (^4].[T7].[TB]). 

Now we can proceed with the stack of iVp 3 = p D3-branes placed at the tip. First 
we require 



at infinity or 



f = k^ -Pifogr* - hge 2 - 1 - j^) . (3.i 



Mt^oo) = Mt^oo) = ^p + Ploge 2 . (3.9) 
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This will fix e as a function of X±, Y^. At the next step we fix Y 5 using (I3.2p with 
^D3 = 0, iVpg = p. Now we only need to fix X\. An easy way to do that is to look 
at the force on probe Nd 3 D3-branes placed at some location on the conifold in the 
presence of N^ 3 D3-branes located at the tip. At first we would like to show that 
the resulting force has the expected form i.e. it is the same at the linear in Nn 3 
and A% 3 order as the force on N D3 probe D3-branes sitting at the tip produced by 
N D3-branes sitting at some location on the conifold and affecting the background 
through the backreaction (13.21) . 

To this end we consider the linearized equation for $± = e 4A ± a where the warp- 
factor h = e~ 4A (here we introduced A = A — p — x/2 to distinguish it from A) and 
a is the RR 4- form C 4 = a dx° A .. A dx 3 . The equations for $ ± are 



V 2 $+ = e" 4A (V$ + ) 2 + e 8A |G+| 2 + e 8A N D3 n <5(location of D3 - branes) , (3.10) 
V 2 $_ = e" 4A (V$_) 2 + e 8A |G_| 2 + e sA N m K ^(location of D3 - branes) . (3.11) 

Here is some constant and we used the notations of [H] for the ISD and 

AISD parts of the 3-form flux G± = (*q ± i)G 3 . In the case of the ISD solution, like 
the KS solution perturbed by some mobile D3-branes, $_ = G_ = and therefore at 
the linear order the equation for the warp-factor is 

V 2 /i = ~\\G+\ 2 - ^iVD35(location of D3 - branes). (3.12) 
Assuming we place D3's at some point Xm we have for the perturbation of h 

8h{X) = ~N m G(X,X m ) , (3.13) 



where G is the Green's function on the conifold [IS]. Now treating the D3-branes as 
probes the corresponding D3-D3 potential at the linear in iV^, A% 3 order is 



Vos-ds = 2T 3 N T53 h-\X 753 ) = KT f { ^ G(X m , X D3 ) . (3.14) 

We can reproduce exactly the same answer (including the coefficient) in linear 
order if we treat the D3-branes as sources changing the geometry. Then it follows 
from (13.1 ip that at linear order 

V 2 $_ = h~ 2 Nv 3 K ^(location of D3 - branes) . (3.15) 
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As a small remark, we emphasize that the mode <£>_ is the only one directly sourced 
by the delta-function coming from the D3-branes [2] . That provides the IR boundary 
condition for ( E35D which follows from fEUOl l3TT|) . 

The solution for the linear mode of is 

kN— 

and the corresponding potential for the probe D3-branes at the linear in Nd 3 , N^ 3 
level is 

Vbs-ua = = ^^G(X D3 , X m ) . (3.17) 

The results ( I3.14p and ( I3.17P are the same because of G(X D3 , Xjj 3 ) = G(Xp 3 , X D3 ). 
Thus we have shown that the force on D3-branes induced by the D3-branes is the 
same as the force on the D3-branes induced by the D3-branes. Clearly, the argument 
above is not based on smearing and is valid in the case of localized branes. 

Now we can use the SU(2) x SU(2) invariant mode of (I3.16P to fix X\. A direct 
calculation gives [H] 

= 2 -e- 2 % . (3.18) 

This should match with the SU (2) x SU (2) invariant mode of (I3.16P with N^ 3 = p 

32Xi ~ 127rP 2 A n o 

$ = Xs = — 3 19) 

- 3(cosh(t) sinh(t) - 1) 2 /3 ' 1 h{0) 2 M 2 ' 1 ; 

This explains why the equation for ^ can be integrated in the form ( I2.26P and fixes 
X\ to be 

12ttP 2 p 

x ' = moFIF • (3 - 20) 

This fixes all other integration constants X i: Yi and hence determines in the unique 
way the lowest KK mode for the linearized solution describing p D3-branes sitting at 
the tip of the deformed conifold. 

The Page charge of the resulting solution is given (mod M) by Q = —^f-pYh = ~P 
in agreement with our expectations. 



1 We thank the authors of [5] who pointed this out to us. 
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4 Tension of D3-branes 



As a check of our solution we will calculate the mass of the metastable state by 
calculating the ADM mass of the corresponding gravity solution found above. Our 
calculation will closely follow a similar calculation carried out in [6]. The necessary 
ingredients are the extrinsic curvature of the eight-dimensional space /C located at 
given time Xq and infinite radius r — > oo and the norm of the time-like Killing vector 
N '. The resulting mass 

E = -\- [ JV7C , (4.1) 
K io J 

can be rewritten in the following form 

M = - 4 ^ - e iA+ * p A' , E= f d 3 xM . (4.2) 



(27r)V 4 

This expression should be calculated at infinite t. It diverges as e 4t//3 ~ r 4 and hence 
should be regularized, presumably by subtracting the tension E of the supersym- 
metric configuration with the same UV asymptotic. Following [T2] we subtract from 
(14.21) the value of the superpotential at infinity 

■M renorm = "TTT^I ( ^ ~ Z) ■ (4-3) 



(2yr) 4 a 

This choice ensures that renorm is zero for any supersymmetric solution that satisfies 
the BPS equations (12.31) i.e. for which all £j vanish. To calculate the mass for the 
configuration in question we expand (14. 3 j) at linear order and notice that only terms 
proportional to survive 

^=(2^«' + «- < 4 - 4 > 

Given that + ^ 4 = A 4 and using (I2.33|) and (I2.36P and also the explicit expression 
for Xi (I3.20p we get for the tension of the meta-stable state at the linear in A% 3 order 

-Mr™ = ^ • (4.5) 

This result is in complete agreement with the probe approximation. Let us emphasize 
here that (14. 5 j) depends on the IR integration constant A4 which is sensitive to the 
IR boundary conditions. Hence our calculation provides a non-trivial check of the 
boundary conditions and the resulting solution. 
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5 Gravity dual for the KS theory perturbed in UV 



As a by product of our analysis we obtained the general solution for the KS back- 
ground perturbed by the SU{2) x SU(2) invariant, scalar, parity and Z-symmetry§| 
even operators W± = W\ ± W$ and with small coefficients. Some of these 

operators are irrelevant and the corresponding modes will destroy the AdS behavior 
at infinity. To avoid this the perturbation at the UV cutoff scale A should be small 
in the conventional units 

where dimensionless e <C 1. 

Since a non-zero X\ indicates presence of the D3-branes we put it to zero. This 
drastically simplifies the resulting solutions. The requirement that the solutions are 
regular in the IR implies X 4 = A 6 = and X s = —PX 5 . At the same time X 5 
can be expressed through X 2 ,X 3 ,X 7 with help of the zero energy condition (I2.29p . 
Moreover we fix the gauge Y\ — to ensure the warp-factor behaves at infinity as in 
the unperturbed case. 

Now we have seven integration constants left. They correspond: X2, X 7 to bottom 
components of W±, Y 8 ,Y 5 to top components of W± i.e. coupling constants; X 3 , Y 7 
and I4 to bottom, medium and top components of WlWl correspondingly. 



5.1 KS with softly broken SUSY 

A particularly interesting example of the UV perturbation discussed above is a small 
mass for gaugino bilinears 2 which softly breaks SUSY. This a relevant perturbation 
of field theory which corresponds to the fluctuation of metric 02 — and the three- 
form flux 5(\(f — k) — F) = |(05 — 6 ) — 07 . The regularity at the IR together with 
the absence of perturbation by the irrelevant (and marginal) operators in the UV 
results in two free parameters X2, V7 

X 8 = -PX 5 = -(2X 2 + 3PX 7 ) , X 6 = X, = X 1 = X 3 = , ^ = 0,1^5.(5.2) 

2 The X-symmetry acts by exchanging the A and B bifundamental fields of the KS theory accom- 
panied by a charge conjugation. See e.g. [19] . 
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We should choose Y 5 = as there is no D3 or D3-branes at the origin. At the same 
time the deformation parameter e should be chosen such that it compensates 05,06 
at infinity. 

Let us point out that a similar solution describing the KS theory perturbed by 
a dimension 3 operator was found in [13]. It corresponds to the one-parametric 
subfamily PXf = —X 2 . 
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